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Abstract. We discuss the effect of structure-preserving perturbations 
on complex or real Hamiltonian eigenproblems and characterize the 
structured worst-case effect perturbations. We derive significant ex- 
pressions for both the structured condition numbers and the worst-case 
effect Hamiltonian perturbations. It is shown that, for purely imaginary 
eigenvalues, the usual unstructured perturbation analysis is sufficient. 



1. Introduction 

For structured eigenvalue problems, algorithms that preserve the underly- 
ing matrix structure may improve the accuracy and efficiency of the eigen- 
value computation and preserve possible eigenvalue symmetries in finite- 
precision arithmetic. The concept of backward stability, i.e., the require- 
ment that the computed eigenvalues are the exact eigenvalues of a slightly 
perturbed matrix, is extended by requiring the perturbed matrix to have the 
same structure as the original one |B87j . In order to assess the strong nu- 
merical backward-stability of the employed structured algorithm and not to 
overestimate the worst-case effect of perturbations, it is appropriate to con- 
sider suitable measures of the sensitivity of the eigenvalues to perturbations 
of the same structure. 

Our analysis is based on the perturbation expansion of a simple eigen- 
value. The structured condition number of an eigenvalue A G A is indeed a 
first-order measure of the worst-case effect on A of perturbations of the same 
structure as A. The structured conditioning measures we deal with can be 
computed endowing the subspace of matrices with the Frobenius norm; see, 
e.g., [TTi3llKKT06llRn6] and references therein. 

Here we are concerned with the Hamiltonian eigenvalue problems. They 
arise from a number of applications, particularly in systems and control 
theory. Algorithms and applications of Hamiltonian eigenproblems are dis- 
cussed, e.g., in [BKM05] and references therein. 

We investigate the sensitivity of the eigenvalues of a complex [real] Hamil- 
tonian matrix with respect to complex [real] Hamiltonian perturbations. 
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Attention has been paid to derive straightforwardly computable formulae, 
for both the structured condition numbers and the maximal structured per- 
turbations, that result to be explicative of the relationship between unstruc- 
tured and structured conditioning. It is well known, see, e.g., |W65j . that the 
worst (unstructured) perturbation which may affect a simple eigenvalue A of 
a given matrix A arises under the action of the matrix yx* , with x and y cor- 
responding right and left eigenvectors, normalized to have \\x\\ 2 = ||y|| 2 = 1- 
As a matter of fact, if we add the Hamiltonian-structure requirement, just as 
for sparsity-structures |NP06] or symmetry-patterns [NP07| . we prove that 
the Hamiltonian matrix which yields the structured worst-case effect per- 
turbation is the structured analogue of yx*, i.e. the normalized projection of 
yx* onto the space of the complex [real] Hamiltonian matrices. Notice that 
a further condition on the angle between the right and left eigenvectors is 
here required in the normalization. 

The paper is organized as follows. Section [2] provides formulae for the 
eigenvalue structured condition numbers and for the maximal Hamiltonian 
perturbations. In Section [3l the real case is addressed. Numerical tests are 
presented in Section SJ Finally, conclusions are drawn in Section 

2. Eigenvalue structured conditioning of Hamiltonian matrices 

We start with some notation and definitions. We denote by Tic the linear 
subspace of 2ra-dimensional complex Hamiltonian matrices, i.e., 

n c = {Q e C 2nx2n : QJ = (QJ)*} 

Q= \L -K*) : K ' L ' M e £nXn ' L = L * ' M = M * 
where J is the fundamental symplectic matrix, i.e., 

/„ 



J 



-h. 



with /„ the n x n identity matrix. We observe that the Hamiltonian matrix 
Q is uniquely determined by the 4n 2 real parameters given by the 2n 2 real 
and imaginary parts of the elements of K, the In (real) diagonal elements 
of L,M, and the 2n(n — 1) real and imaginary parts of the elements of the 
strict lower triangular parts of L, M. 



Proposition 2.1. The closest Hamiltonian matrix to a given matrix A 
A x A3 
A 2 A A/ 



A\ A'i \ 

' G (£2nx2n yjjjfa reS p ec f f the Frobenius norm is 



If A x -A* d A 3 + A*,\ 1,. T .* T s 
A ^ = 2{a 2 + 4 a\-4) = 2 {A + JAJ) - 

Proof. For Q = ^ ^ we have 

\\A - Qf F = \\K - A-xfp + \\K + Alf F + \\L - A 2 \\ 2 F + ||M - A 3 f F 
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Since L, M are only requested to be Hermitian, it is known that the last 
two terms in the right-hand side are minimized for L = ^(A 2 + ^2) an d 
M = I (.A3 + A3). On the other hand, the minimization of 

G{K) = \\K-A 1 f F + \\K + Al\\ 2 F 

leads to K = \{A\ — A|). The proposition is thus proved. □ 

To state the next results we need some additional notation and defini- 
tions. It will be useful the following characterization of the space of complex 
Hamiltonian matrices, 

H c = {H + iW: H G H, W€W}, 

where H [resp. W] is the linear space of real Hamiltonian [resp. skew- 
Hamiltonian] matrices, i.e., 

U = {H£ M 2 ™x 2 ™: HJ = {HJ) T } , 

W = {W G R 2nx2n : WJ = -{WJf} . 

In particular, if A = ^R(A) + r3(A) then A\ Hc = $t(A)\ n + i9(A)| w where, 
for B G R 2 ™ x2 ™, 

B\ H :=±(B + JB T J) B\ W :=^(B-JB T J). 
We finally introduce the normalized projection of a matrix A onto He as 

\\ a \h c \\f 

Theorem 2.2. Let X be a simple eigenvalue of a Hamiltonian matrix Q, 
with corresponding right and left eigenvectors x and y normalized to have 

||z|| 2 = ||y|| 2 = l, %(y*Jx) = 0. (2.2) 

Given any Hamiltonian matrix E with \\E\\p = 1, let Ac(i) be an eigenvalue 
of Q + tE converging to X as t — > 0. Then, 



|A E (0)| < max 
Moreover, 



y*Gx 



y*x 



\G\\ F = i,GeH C( = m 't If ■ (2.3) 

\y*x\ 



\ E (<S)= m fan: if E = ± yx *u. (2.4) 

y*x 

Remark 2.3. By standard perturbation theory, the variational estimate in 
(|2.3p holds for any choice of the right and left eigenvectors x and y. Instead, 
the assumption (|2.2p is needed for the explicit expressions of the maximum 
and of the maximizers to hold. 

Proof of Theorem \2.2[ We shall use the real scalar product (A, B) := 
Trace(A T 5), A,B e R 2nx2n , note \\A\\ F = y/{A,A). Our target is to 
compute 

HTgmax{\y*(H + iW)x\: H eH, W G W , \\H\\ 2 F + \\W\\ 2 F = l} . 
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The requirement H £ % means that H J is symmetric, i.e., it is orthogonal 
(w.r.t. the above scalar product) to the space of skew-symmetric matrices. 
Similarly, W £ W means that W J is orthogonal to the space of symmet- 
ric matrices. Therefore, we study the auxiliary (unconstrained) variational 
problem for the function 

G(H,W) = \y*(H + iW)x\ 2 + X(\\H\\ 2 F + \\W\\ 2 F -l) + (A,HJ) + (T,WJ) , 

which depends on 4n 2 + 1 Lagrange multipliers, precisely the real parameter 
A, the skew-symmetric matrix A, and the symmetric matrix V. Such param- 
eters will be eventually determined by requiring that the solution satisfies 
all the constraints. 

It is now useful to introduce the matrices £ = 5i(yx*), r\ = Q(yx*), so 
that 

yx* = £ + i 77 , yx T = £ — i 77 

and therefore 

y*(H + iW)x = Trace({yx T ) T (H + iW)) = Trace((£ T - ir, T )(H + i W)) 
= (£,H) + ( V ,W) +i((£,W) - ( V ,H)) . 

Whence, 

G(H, W) = «£, H) + ( V , W)f + «£, W) - (v, H)) 2 

+ \{\\H\\ 2 F + \\W\\ 2 F - 1) + (A, HJ) + (T, WJ) , 

and, by (I2.2p . also noticing that y* Jx = Trace(xy* J) = (£, J) — i (77, J), 

||£||! + ||r,||! = l, ( V ,J) = 0. (2.5) 

We now look for the critical points of G. The gradients V hG(H,W), 
X7wG(H, W), expressed in the form of 2n x 2n matrices, read, 

V H G(H, W) = 2ai - 2br] + 2XH + AJ T , 

V W G(H, W) = 2ar] + 26£ + 2XW + TJ T , 

where we introduced the scalar functions, 

a = {£,H) + {ri,W) , 6 = (£, W) - ( V , H) . 

The critical points (H, W) are therefore solutions to 

a /- b 1 A tT 

A A 2A , 2 _gx 

By imposing HJ = {HJ) T and VKJ = -(WJ) T , using that A = -A T , 
T = T T , and recalling that JJ T = I, J T = — J, we get, 

A = -a(£J + J£ T ) + b(T]J + Jri T ) , r = —a(rjJ - Jrj T ) - b(£J - J£ T ) . 



STRUCTURED CONDITIONING OF HAMILTONIAN EIGENVALUE PROBLEMS 5 

Plugging these expressions in (12. 6ft and calling a = —X~ 1 a, f3 = X~ x b, we 
conclude that the critical points have to be of the following form, 

H = a{i + J? J) + P(n + Jr, T J) , 

W = a(<n- Ji] T J) - /?(£ - J^J) ■ 

To determine the (still) unknown real parameters (a,/3), it remains to max- 
imize \y*(H + iW)x\ 2 with the constraint + ||W|||, = 1 for (H,W) as 
in (|2.7jl . To compute these quantities we first notice that, as (£ T ± ir/ T ) J 
are rank-one (complex) matrices, 

Trace((^ T + ir ? T )J(^ T ±ir ? r )J) = Trace((£ T + irj T ) J)Trace(£ T ± ir/ T ) J) , 
from which, by identifying the real and imaginary parts, we get 

Trace (£ T J£ T J) = TYace(£ T J) 2 , Trace (r] T Jrj T J) = Trace (r? T J) 2 , 
Trace (£ T Jr/ T J) = Trace (r? T J£ T J) = Trace (r? T J) Trace (£ T J) . 



Therefore, by (1231) . 

(e, ^ T ^> = <e, </> 2 , fa, V^} = (v, j) 2 = o, 
<e, V j> = fa, = <e, fa, j) = o , 

so that, one easily computes, 

\y*(H + iW)x\ 2 = (1 + (£, J) 2 ) V + (1 - <£, J) 2 ) 2 /3 2 , 

II^IIf + II^IIf = 2(1 + <£, J> V + 2(1 - <£, J) 2 )/? 2 . (2.8) 

Moreover, as yx*\% [resp. (iyx)*\%\ is given by (|2.7|) for (a,/3) = (1/2,0) 
[resp. (a,/3) = (0, 1/2)], by flZE) we get 

, 2 _ i + (e,j) 2 w I|2 _ i-<e,j> 2 



I^1wIIf = ^ — . II ^y x )*\n\\ F 



2 ' llv " ' ' L " r 2 

If ||(iyx)*|^|| > 0, i.e., | (£, J) | < 1, the (constrained) maximum of \y*(H + 
i V^):r| 2 is reached for a 2 = [2(l+(£, J) 2 )]" 1 and ft 2 = 0, so that the maximal 
perturbation is obtained for 

-E± = ± ; ± l - = = ±yx \m . 

2(i + (e,j) 2 ) V 2 ( 1 + ^ J ) 2 ) 

On the other hand, if ||(i?/x)*| w || = then a 2 = [2(1 + (£, J) 2 )] -1 and 
H + i W = E± holds trivially. In both cases, 



\y E±x\ = y = Hs/sc • 

The theorem is thus proved. □ 
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Remark 2.4. By Theorem 12.21 the structured condition number is explicitly 
given by 



which can be also used for the structured backward error analysis of the 
algorithms in [BMX98, BK05]. Since the (unstructured) condition number 
of a simple eigenvalue is k(A) = we have k-^ c (A) = k(X) if 

yx* is Hamiltonian. As shown in the next proposition, this can occur only 
if A is a purely imaginary eigenvalue. 

Proposition 2.5. Let X be a simple eigenvalue of a Hamiltonian matrix Q 
with corresponding right and left eigenvectors x and y normalized such that 
\\ x \\2 = IMk = 1- If yx* is a Hamiltonian matrix or a skew-Hamiltonian 
matrix (i.e., {yx*J) = —(yx*J)*) then 3?A = 0. 

Proof. We prove the statement in the case yx* is a Hamiltonian matrix, the 
case yx* is skew-Hamiltonian can be treated similarly. We first show that 
if yx* is Hamiltonian then y = ±Jx. Indeed, since J T = — J, the identity 
yx* J = (yx* J)* can be recast as yk(Jx)j = (Jx)ky~j for any k,j = l,..., 2n. 
As y*x ^ 0, there exists at least one index k for which (Jx)k,yk 7^ 0. Then, 
it is readily seen that the above identity implies Jx = cy for some c £ R. 
As \\x\\2 = H2/H2 = 1; we conclude that y = ±Jx. It remains to show 
that y = ±Jx implies 3ftA = 0. To this end, it is enough to notice that, 
as Q is Hamiltonian and J 2 = —1, the identity Qx = Xx can be recast 
as (Jx)*Q = —X(Jx)*. Since also y = ±Jx is solution to y*Q = Xy* , we 
conclude that A + A = 0. □ 

3. Structured conditioning of real Hamiltonian eigenproblems 

In several situations, we are interested in real structured perturbations of 
real Hamiltonian matrices. For these purposes, we are led to consider the real 
structured condition number ac%(A), where only real Hamiltonian perturba- 
tions are taken into account. For any smooth manifold S, an explicit for- 
mula for the real structured condition number ks(X) is given in |KKT06l Eq. 
(2.16)]. Unfortunately, such formula, as its counterpart in the complex case, 
requires the characterization of the tangent space of S and the construction 
of a pattern matrix relevant to S. On the other hand, it yields the following 
bound, 



Actually, we can give a more explicit expression of k-^(A), which is the 
content of the following theorem. 

Theorem 3.1. Let X be a simple eigenvalue of a real Hamiltonian matrix 
Q with corresponding right and left eigenvectors x and y normalized to have 



x\\ 2 = \\y\\ 2 = 1 , %(y*Jx) ^(y*Jx) = Trace(K(ya;*) J $s(yx*)) . (3.1) 




V2 



< K H {X) < KH C W ■ 
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Given any real Hamiltonian matrix E with \\E\\p = 1, let A^(i) be an eigen- 
value of Q + tE converging to A as t — > 0. Then, 



|A#(0)| < max 



y*Hx 



max 



# f = l,Hen 

y^x 

{ \Myx*)\ n \\ F _ \myx*)\ n \\F} 



V 



|2/ x| |y x| 



(3.2) 



Moreover: 

j; i/||K(7/x*)| w || F > ||3?(yx*)W| F iften 



t/*x l|K(yx*)| w )|| 
S; ||fc(l/z*)|w||jr < \\%(yx*)\ H \\F then 

X E (0) = ^MZ for E - . 



f 



5; i/||»(i/s*)|«||F = ||9(i/x*)|«||F ^en 

t / n x - „ ^(7/x*)| w cos g + 5(j/x*)|^ sing 

Af(0 = for any E = — — -— , 

y*x W(vx*)\h\\f 

with 6 G [0, 2vr]. 

Proof. Our target is to compute 

argmax [\y*Hx\ : H EH , \\H\\ F = l} . 

As before, we introduce the real matrices £ = 9?(yx*), ?? = Q(yx*). By 
arguing as in the proof of Theorem 12.21 we find that the argmax have to be 
of the following form, 

H = a(£ + J£ T J) + /0(7j + Jn T J) . (3.3) 

As (f + J£ T J, 77 + Jr] T j) = 2 (£, 77) + 2 (£, J) (77, J) and recalling y* Jx = 
(£, J) - i (r/, J), by fl^T]) we have, 

lieill + IMIf = 1 , (£ + ^ V + Jr, T J) = . 

Whence, a straightforward computation shows that if H = a(£ + J^ T J) + 
+ J if J) then 

|t/*#x| 2 = £>fa 2 + L> 2 /3 2 , ||iT||| = 2,05a 2 + 2D^ 2 , 

where 

1 



d* ■■=^M+jej\\ F = M\\F+(z,j) 2 
1 



A, := §lk + ^ T J||| = \\rj\\ F + (n, J} 2 
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Moreover, as $l(yx*)\u [resp. ^s(yx*\%)] is given by (|3.3p for (a, /3) = (1/2,0) 
[resp. (a, 13) = (0, 1/2)], it follows that 



Myx*)\H\\ 2 F = ^ 



myx*)\n\\ 2 F = % 



If both and D v are positive the (constrained) maximum of \y*Hx\ 2 
is reached for (a 2 ,/3 2 ) = ((2 J D C )~ 1 ,0) if £> 5 > D„ and for (a 2 ,/? 2 ) = 
(0, (2D V )- V ) if < D n . In the critical case -Dg = -D^ =: D, the func- 
tion |y*i?x| 2 is constant and any pair (a 2 ,/3 2 ) satisfying the constraint 
2D[c? + /3 2 ) = 1 is allowed. Therefore, the maximal perturbation is ob- 
tained for 



E v = ± 



i] + Jrj T J 



if > D v , 



if D^<D V , 



while, in the critical case -Dg = =■ D, for 

(f + J£, T J) cos + (77 + Jrj T J) sin 



'2D 



G [0, vr] 



On the other hand, if = [resp. = 0] then a 2 = (2D^) 1 [resp. 
(3 2 = (2D V )~ 1 ] and H = [resp. if = £,,] holds trivially. In all cases, 



'Dt 



\y*E v x\ 
The theorem is thus proved. 



□ 



Remark 3.2. We observe that yx* is Hamiltonian if and only if Q(yx*)\y_ = 
and K(yx*) = $t(yx*)\u. Likewise, yx* is skew-Hamiltonian if and only if 
^(u x *)\h = and Q(yx*) = $s(yx*)\u. Therefore, by Theorem 13.11 



(\myx*)h 



I y x 
\Q(yx* 

\y*x\ 



if yx* is Hamiltonian , 

if yx* is skew-Hamiltonian . 



By |BK031 Theorem 2.1] it follows that Ky,(X) is equal to the real (unstruc- 
tured) condition number. Moreover, by Proposition 12.51 this occurs only if 
A is a purely imaginary eigenvalue. 

Remark 3.3. If is worthwhile to remark that if the eigenvalue A is real, since 
the right and left eigenvectors are also real, then (|2.2I) and (|3.ip are trivially 
satisfied and ku{\) = nu c (A). 
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0.3128 



0.3127 



0.3126 



0.3125 



0.3124 




0.3203 



0.3204 



0.3205 



0.3206 



0.3207 



0.3208 



Figure 1. Example 4.1. Perturbations induced in the eigen- 
value A by the complex structured maximal perturbations 
(green *), the real structured maximal perturbations (red — ) 
and the unstructured maximal perturbations (blu — ). 



4. Numerical tests 

Example 4.1. Consider the Hamiltonian matrix H 
where 

'0 2\ _ (-0.1 o.r 

0/' 2 ~ I 0.1 0.1 



Hi H3 
H2 H4 





0.1 



0.1 
-0.1 



»4x4 



It is well known that if A € C — (RUiR) is an eigenvalue of a real Hamiltonian 
matrix, then also —A, A, —A are eigenvalues, since the spectrum is symmetric 
with respect to the real and the imaginary axis. According to Theorems 12.21 
and l3.ll one derives for the eigenvalues of H the condition numbers reported 
in Table[TJ Let us focus on the eigenvalue A = 0.3205+i 0.3126, marked by a 



A 


«(A) 






±0.3205 ±i 0.3126 


2.3513 


1.6631 


1.5851 



Table 1. Example 4.1. Unstructured, structured and real 
structured eigenvalue condition numbers. 
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black point in Figure[U and consider e-norm perturbations, i.e. perturbation 
matrices of Frobenius norm equal to s, with e = 10 -4 . 

The complex structured maximal perturbation matrices in (I2.4p shift A 
to one of the points marked by a green star in Figure H) As a matter of 
fact, with respect to structured perturbations, the eigenvalue is not equally- 
sensitive in any direction of the complex plane, since its eigenvectors have 
been normalized such that (I2.2j) is satisfied. 

As for the real structured perturbations, holding in this case the equality 
ll^(y a; *)lwl|.F = ll^y^Olwll-F) the maximal perturbation matrices are given 
by case 3) in Theorem 13.11 and they shift A to the boundary of the red circle 
of radius K%{X)e. 

We point out that the unstructured worst-case effect perturbation matri- 
ces ee ld yx*, for 9 € [0, 2tt], shift in their turn A to the boundary of the blue 
circle of radius n(\)e. Analogous results hold for the other three eigenvalues 
of H. 

Example 4.2. Consider the same Hamiltonian test matrix as in Example 
4.1, except for the block H3 that is here multiplied by the factor 0.1. We 



A 


k(X) 






±0.1786 ±i 0.1771 


7.0263 


6.3769 


6.3184 



Table 2. Example 4.2. Unstructured, structured and real 
structured eigenvalue condition numbers. 



consider e-norm perturbations in the matrix and analyze the effects induced 
in A = 0.1786+i 0.1771. The complex [real] structured maximal perturbation 
matrices shift A to one of the points marked by a green star [red circlet] in 
Figure[2j Since < ||9f(ya;*)|^||j?, case 2) in Theorem 13.11 holds ; 

nevertheless, for the sake of completeness, let us take into account even the 
(non-maximal) real Hamiltonian matrices 

= ?R.(yx*)\ncos9 + $s(yx*)\ n sm6 . . 

6 \\&(yx*)\ n coaO + Ql(yx*)\ H 8m0\\ F ' { ' } 

and notice that the perturbation matrices e Eg, for 9 € [0, 2tt], have the effect 
of shifting A to the boundary of the ellipse in Figure [H whose major axis 
gives indeed the direction induced by the worst-case effect perturbations. 

5. Conclusions 

We have derived simple expressions for structured eigenvalue condition 
number and for the worst-case effect perturbation in case of complex and 
real Hamiltonian matrices, which turns out to be a rank-2 complex matrix 
in the former and a rank-4 real matrix in the latter. In Proposition 12.51 and 
Remark 13.21 we showed that the structured and unstructured eigenvalue 
condition numbers are equal in case of purely imaginary eigenvalues, the 
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0.1778 
0.1776 
0.1774 
0.1772 
0.177 
0.1768 
0.1766 



0.1778 0.178 0.1782 0.1784 0.1786 0.1788 0.179 0.1792 0.1794 



Figure 2. Example 4.2. Perturbations induced in A by the 
complex structured (green *), real structured (red o) and 
unstructured (blu — ) maximal perturbations, and by the real 
structured perturbations in (|4.ip (red — ). 

unstructured worst-case effect perturbations resulting to be Hamiltonian. 
On the other hand, numerical tests show that the structured perturbation 
analysis may be of some interest for eigenvalues not lying in the imaginary 
axis; see, e.g., the tests reported in Section |H 
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